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Abstract 

We present a simple scheme for taking into account the resonant continuum coupling in the 
Relativistic Mean Field- BCS (RMF-BCS) calculations. In this scheme, applied before in non- 
relativistic calculations, the effect of the resonant continuum on pairing correlations is introduced 
through the scattering wave functions located in the region of the resonant states. These states are 
found by solving the relativistic mean field equations with scattering-type boundary conditions for 
continuum spectrum. The calculations are done for the neutron-rich Zr isotopes. It is shown that 
the sudden increase of the neutron radii close to the neutron drip line, the so-called giant halo, is 
determined by a few resonant states close to the continuum threshold. 
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I. INTRODUCTION 

As recognized long time ago the basic features of the superfluidity are the same in 
atomic nuclei and infinite Fermi systems. Yet, in atomic nuclei the pairing correlations have 
many features related to the finite size of the system. The way how the finite size affects 
the pairing correlations depends on the position of the chemical potential. If the chemical 
potential is deeply bound, like in stable and heavy nuclei, the finite size influences the pairing 
correlations mainly through the shell structure induced by the spin-orbit interaction. The 
situation becomes more complex in nuclei close to the drip lines, where the chemical potential 
approaches the continuum threshold. In this case the inhomogeneity of the pairing field 
produces strong mixing between the bound and the continuum part of the single-particle 
spectrum. Due to this mixing the quasi-particle spectrum becomes dominated by resonant 
quasi-particle states, which originate both from single-particle resonances and deep hole 
states J 1Q- 

The continuum effects on pairing correlations is commonly taken into account in the 
Hartree-Fock-Bogoliubov (HFB) (| or Relativistic-Hartree-Bogoliubov (RHB) []| approach. 
In most of these calculations the continuum is replaced by a set of positive energy states 
determined by solving the HFB or RHB equations in coordinate space and with box bound- 
ary conditions 3,0]. Due to this fact the genuine continuum effects, as the widths of the 
quasi-particle states, are not accounted for in these type of calculations. 

Recently the HFB equations were also solved with exact boundary conditions for the 
continuum spectrum, both for a zero range and a finite range pairing force [10]. It 
was thus shown that close to the drip lines the discretisation of the continuum generally 
overestimate the pairing correlations. A similar conclusion was obtained earlier in a simpler 
BCS approach, in which the resonant part of the continuum was studied Comparing 
these BCS results with the exact HFB solutions P] one finds that the effect of the continuum 
on pairing correlations is given mainly by a few resonant states close to the continuum 
threshold. 

For the relativistic models an exact solution of the continuum spectrum is not available 
yet, neither for RHB nor for RMF-BCS approach. A comparison between RHB and RMF- 
BCS calculations, performed by using box boundary conditions, is discussed in Ref. 
This comparison indicates also the special role played by the resonant states, which in these 
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calculations are approximated by positive energy states. This approximation works well 
only if the positive energy states correspond to very narrow resonances. Moreover, since a 
discrete representation of the continuum does not provide a direct measure of the width of 
the resonant states, the selection of the relevant positive energy states is ambiguous if the 
resonant states close to the continuum threshold are not very narrow. 

The scope of this paper is to show how the resonant continuum can be treated accurately 
in the RMF-BCS approach. The single-particle states belonging to the resonant part of the 
continuum spectrum will be calculated by solving the RMF equations with scattering-type 
boundary conditions. Then the resonant continuum will be handled in the BCS equations in 
a similar way as in the non-relativistic HF-BCS calculations [12]. This approach is applied 
for the case of Zr isotopes, for which earlier calculations predict a very large neutron skin 
close to the neutron drip line. It is shown that the sudden increase of the nuclear radii in 
these isotopes is essentially determined by a few single-particle resonant states close to the 
continuum threshold. 

The article is organized as follows. In Section 2 we discuss shortly the scattering type so- 
utions of the relativistic mean field equations and we introduce the resonant-BCS equations 



12j. Then in Section 3 we present the results of the calculations for Zr isotopes. Section 4 



contains the summary of the paper. 

II. RESONANT STATES IN THE RMF-BCS APPROACH 
A. Continuum-RMF solutions 

In the relativistic mean field approach the nuclear interaction is usually described by the 
exchange of three mesons: the scalar meson a, which mediates the medium-range attraction 
between the nucleons, the vector meson lu, which mediates the short-range repulsion, and 
the isovector-vector meson p, which provides the isospin dependence of the nuclear force. 



The equations of motion are commonly derived from the effective Lagrangian density 
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-\f, v F^ - e^il^U'V , (1) 

where a non-linear self-coupling is considered both for a and uj mesons. The vector fields H, 
G and F are given by 

H^y = d^ v - d v u^ 
F\iv = d fjL A u — d v A^ 

The nucleons are described by the Dirac spinor field ip, which in the case of spherical 
symmetry can be written as: 

r \Fa-r y jWi ' 



where [V?7m denotes the spinor spherical harmonics while G and F are the radial wave func- 
tions for the upper and lower components, respectively. They satisfy the radial equations: 

^y- + -G — {M + E + V S — V V )F = (3) 
ar r 

- -j- + -F + (M - E + V S + V V )G = (4) 
ar r 

where V s and V v are the scalar and the vector mean fields and k is given by: 



K 



+ if J = / + l/2 
+1 if j = 1-1/2 



(5) 



At large distances, where both the scalar and the vector mean fields are zero, the radial 
equations can be written in the form: 

^ + (a *_^ )G .0 (6, 

where a 2 = E 2 — M 2 . These equations are suited for fixing the scattering- type boundary 
conditions for the continuum spectrum. They are given by: 

G = Car[cos(5)ji(ar) — sin(S)ni(ar)} (8) 



Ca 2 r 

[cos(8)ji„i(ar) - sin(5)ni^ 1 (ar)], (9) 



E + M 

where ji and n; are the Bessel and Neumann functions and S is the phase shift associated 
to the relativistic mean field. The constant C is fixed by the normalisation condition of the 
scattering wave functions and the phase shift 5 is calculated from the matching conditions. 
In the vicinity of an isolated resonance the derivative of the phase shift has a Breit-Breit 
form, i.e., 

dS(E) _ r/2 
dE ~ (E r - E) 2 + T 2 /4 1 1 

from which one estimates the energy and the width of the resonance. In the vicinity of a 

resonance the radial wave functions of the scattering states have a large localisation inside 

the nucleus. Close to a resonance the energy dependence of both components of the Dirac 

wave functions can be factorized approximatively by a unique energy dependent function 

Q|. As in the non-relativistic case 3|> this energy dependent factor is the square root of 

the Breit-Wigner function written above, or, equivalently, the square root of the derivative 

of the phase shift. Using this property all the matrix elements of a two-body interaction 

between these scattering states can be expressed in term of a unique matrix element, i.e. the 

one corresponding to the scattering state with energy equal to the energy of the resonance. 

This property is employed below for the treatment of the resonant continuum in the BCS 

equations. 

B. Resonant states in the BCS approach 

Since the meson exchange forces are not describing properly the pairing correlations in 
nuclei, the relativistic mean field is combined usually with non-relativistic pairing models. 
Here we use for the pairing treatment the BCS approach. The extension of the BCS equations 
for taking into account the resonant continuum was proposed in Refs. [ill . I12I ] . For the case 
of a general pairing interaction these equations, referred below as the resonant-BCS (rBCS) 
equations, reads |l2j |: 

A i = V n.u"j r j + £ V H,ue^ / 9u{e)u u {e)v u {e)de , (11) 
A, = £ V^vx&Fi + E Ke^ye^ I gAe'W(e')vAe')de' , (12) 
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E^ + E/ 9l{e)v 2 v {e)de . (13) 



Here A, are the gaps for the bound states and A u are the averaged gaps for the resonant 
states. The quantity gl{e) = 2j '^ +1 ^ is the total level density and 5 V is the phase shift of 
angular momentum {l v j v )- The factor g„{e) takes into account the variation of the localisa- 
tion of scattering states in the energy region of a resonance ( i.e., the width effect) and goes 
to a delta function in the limit of a very narrow width. The interaction matrix elements are 
calculated with the scattering wave functions at resonance energies and normalised inside 
the volume where the pairing interaction is active. For more details see Ref . jl^ . 

The rBCS equations written above are applied here with the single-particle spectrum of 
the RMF equations. For the pairing interaction we use in the next section a delta force, i.e. 
V = V()5(fi — r 2 ). In this case the matrix elements of the pairing interaction are given by: 

[nn) 0+ \V\ (r 2 r 2 ) 0+) = ^ I dr± (<% G T2 + F T * F T2 f (14) 

For the resonant states these matrix elements are calculated as mentioned above, i.e. using 
the radial wave functions evaluated at resonance energies and normalised inside a finite 
volume. 

The RMF and the rBCS equations are solved iteratively. At each iteration the densities 
are modified through the occupation probabilities provided by the rBCS, as in the non- 



relativistic HF-rBCS calculations 



121. 



III. RMF- RBCS CALCULATIONS FOR NEUTRON-RICH ZR ISOTOPES 

Zr isotopes were discussed recently in connection to the so-called giant halo structure, 
which these isotopes may develop close to the neutron drip line ^|. In Ref. 17 1 these 
isotopes were calculated by solving the RHB equations in coordinate representation and 
using box boundary conditions. The mean field was described by the parameter set NLSH 
jisl ] and for the pairing interaction was employed a density dependent delta interaction. In 
the calculations were introduced all the positive energy states up to 120 MeV. 

In order to compare our calculations with the RHB predictions of Ref.Q] we use for 
the mean field the same parameter set, i.e. NLSH. The results are not much different even 



when we use other parameter sets, e.g. NL3 and TM1 



The appropriate choice for 
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the pairing interaction is more difficult because the pairing correlations estimated with a 
zero range force depend strongly on the energy cut-off, which is very different in the two 
calculations. Thus in the RMF-rBCS approach we include from all the continuum only 
a few resonant states close to zero energy while in the RHB calculations the pairs are 
virtually scattered in all the positive energy states up to the energy cut off, i.e. 120 MeV. 
This energy cut off is here much larger than the maximum quasi-particle energy calculated 
in RMF-rBCS, which corresponds to the single-particle bound state lsi/2- Due to these 
facts we cannot compare meaningfully the results of the two calculations if we use the same 
zero range force. The best we can do is to choose in the RMF-rBCS calculations a pairing 
force which provides in average pairing energies close to the RHB values, at least for some 
isotopes. Following this procedure we chose in the RMF-rBCS calculations a zero range 
pairing force given by V = Vo(ri — r 2 ), with Vq—- 275 MeV. We found that the strength of 
the force can be actually increased up to about Vo=-350 MeV without affecting significantly 
the separation energies and the nuclear radii shown below. 
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FIG. 1: The energies of bound and resonant single-particle states close to the continuum threshold 
in Zr isotopes. The Fermi energy is shown by the dashed line. 
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First we analyse the behaviour of the single-particle states in the vicinity of the 
continuum threshold. These states play the major role in the formation of the neutron skin 
structure discussed below. In Figure 1 are shown the neutron single-particle levels for the 
Zr isotopes closest to the neutron drip line, i.e. from the mass number A = 120 up to 
A = 138. The dashed line represents the chemical potential, which stays close to zero from 
A = 124 to A = 138. The positive energies shown here are the energies of the resonant 
states. As seen in Figure 1, the states 2/ 5 / 2 , l/ig/2 and H13/2 remain resonant states for all 
the isotopes and their energies are changing with the neutron number in the same way as 
the energy of the bound state l/in/2- The other three states, i.e. 2/7/2, 3p3/ 2 and 3pi/ 2 , are 
resonant states for A < 126 and become loosely bound states for heavier isotopes. Their 
radial wave functions for A = 124 are shown in Figure 2. The figure shows the upper 
components of the radial wave functions calculated at the resonance energies, for which 
their localisation inside the nucleus is the largest. 
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FIG. 2: The radial wave functions of the resonant states 2/ 7 / 2 , 3p 3 / 2 and 3j?i/ 2 in 124 Zr. The plot 
represents the upper components of the radial wave functions calculated at the resonance energies. 
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FIG. 3: The width of the single-particle resonant states shown in Figure 1. 

The widths of the resonant states are plotted in Figure 3. The resonant states lhg/ 2 and 
2/5/2 have rather small widths for all the isotopes. On the other hand the widths of the 
resonant states 3pi/ 2 and 3p3/ 2 are changing dramatically with the neutron number. This is 
especially the case for the resonant state 3pi/ 2 . Considering that these states with low (lj) 
values have a major contribution to the formation of the neutron skin, their wave functions 
should be calculated accurately, both for positive and negative energies. The resonant states 
with high (lj) values give instead the dominant contribution to the pairing correlations. Since 
these resonant states have a small width they could be eventually treated like quasi-bound 
states in the pairing calculations. 

Next we examine the two- neutron separation energies S 2n , i-e., 

S 2n (Z, N) = B(Z, N) - B(Z, N — 2) (15) 



Their values are shown in Figure 4. The empirical values correspond to Ref. 



while the 



RHB results are the one given in Ref.Q. One can see that RMF-rBCS gives practically 
the same results as the RHB calulaculations. The two-neutron separation energies remain 



close to zero all the way from A=124 to A=138, which in RMF would corresponds to the 
filling of the group of states 2/7/2, 3p 3 / 2 and 3pi/ 2 . 
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FIG. 4: Two- neutron separation energies of even Zr isotopes as a function of the mass number A. 

In order to see the amount of the pairing correlations in these isotopes we plotted in Figure 
5 the pairing correlation energies, i.e. the binding energies referred to the RMF values. The 
pairing correlation energy curve shows a minimum for A = 136, which corresponds to the 
filling of the states 2/7/2, 3ps/2- These states have almost the same energy and behave like 
a closed major shell for the pairing correlation energy. 

As discussed in Ref. (5, Q], the pairing correlations become usually stronger when the 
widths of resonant states is not taken into account, i.e. when the resonant states are con- 
sidered as quasi bound states. This effect can be also seen in Figure 5, even if in the present 
case the differences in correlation energies are rather small. 

The most interesting phenomenon in these isotopes is the behaviour of the neutron radii, 
which are shown in Figure 6. First one notices that the RMF-rBCS results follow again 
very closely the RHB values. As shown in Figure 6, the neutron radii increase sharply from 
A = 122 to A = 124. This is mainly due to the filling of the states 3p3/ 2 and 2/7/2, which 
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in A = 124 are very narrow resonances close to the continuum threshold. For A = 124 
the resonant state 3pi/2 has almost the same energy as the resonant state 3p^/2 but its 
contribution to the radius is reduced because it has a rather large width. 




FIG. 5: Pairing correlation energies of even Zr isotopes as a function of the mass number A. 

The behaviour of the nuclear radii close to the drip line is very sensitive to the relative 
occupancy of the loosely bound states and the low-lying narrow resonances with high angular 
momenta. Thus if the diffusivity of the Fermi sea is increasing the pairs scatter from the 
loosely bound states to the narrow resonances, which are more localized around the nucleus. 
Consequently the nuclear radii might decrease if the average pairing gap is increasing. 

This effect can be also seen in the present calculations. In the RHB calculations the 
occupancy of the narrow resonances with high angular momenta is larger than in RMF-rBCS 
calculations. Accordingly, as seen in Figure 6, the RHB radii are smaller than in RMF-BCS 
calculations. This situation is quite general since in the RHB or HFB calculations, based on 
a big energy cut off, the Fermi sea is usually more diffuse than in rBCS-type calculations, 
even if the pairing correlation energies might be rather similar in the two calculations. 
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FIG. 6: Root mean square neutron radii of even Zr isotopes as a function of the mass number A. 
IV. CONCLUSIONS 

In this paper we discussed how the resonant states can be treated accurately in the 
RMF-BCS approach. The resonant states are described through the scattering states lo- 
cated in the vicinity of the resonance energies. These states are calculated by solving the 
RMF equations with scattering-type boundary conditions for the continuum spectrum. In 
the RMF-BCS equations the matrix elements of the pairing interaction involving resonant 
states are calculated by using the scattering states evaluated at the resonance energies and 
normalised inside a finite region close to the nucleus. The variation of the matrix elements of 
the pairing interaction due to the widths of the resonant states is taken into account by the 
derivative of the phase shift. This approximation scheme, used previously in non-relativistic 
HF-BCS calculations, is applied here for the neutron-rich Zr isotopes. It is shown that the 
sudden increase of the neutron radii close to the neutron drip line depends on a few res- 
onant states close to the continuum threshold. Including into the RMF-BCS calculations 
only these resonant states one gets for the neutron radii and neutron separation energies 
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practically the same results as in the more involved RHB calculations. 
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